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Abstract. We consider a system of ordered cars moving in K from right to left. Each 
car is represented by a point in M; two or more cars can occupy the same point but 
cannot overpass. Cars have two possible velocities: either or 1. An unblocked car 
needs an exponential random time of mean 1 to pass from speed to speed 1 (slow- 
to-start). Car i, say, travels at speed 1 until it (possibly) hits the stopped car i — 1 
to its left. After the departure of car i — 1, car i waits an exponential random time 
to change its speed to 1, travels at this speed until it hits again stopped car i — 1 and 
so on. Initially cars are distributed in E according to a Poisson process of parameter 
A < 1. We show that every car will be stopped only a finite number of times and that 
the final relative car positions is again a Poisson process with parameter A. To do that, 
we relate the trajectories of the cars to a M/M/l stationary queue as follows. Space 
in the traffic model is time for the queue. The initial positions of the cars coincide 
with the arrival process of the queue and the final relative car positions match the 
departure process of the queue. 



PACS numbers: 1315, 9440T 



A slow-to-start traffic model related to a M/M/l queue 



2 



1. Introduction 

Interest in traffic models is old. In 1935 Greenshield (see Helbing [14] and Chowdhury et 
al j9] for a review on traffic models) introduced the first traffic study. In 1959 Greenberg 
[13J called the attention to the importance of the area. In 1992 Nagel-Schreckenberg [17] 
study traffic using probabilistic cellular automata, computer simulations and mean field 
models. Slow-to-start models intend to capture the behavior of cars that come out of a 
traffic jam: a driver needs a moment to speed the car. Nagel and Schreckenberg [17], 
Schadschneider and Schreckenberg [19], Gray and Griffeath [12], Chowdhury, Santen 
and Schadschneider [9] and Yaguchi [21] studied these type of models. Other studies of 
traffic have been done by Krug and Spohn [15], Barlovic, Santen, Schadschneider and 
Schreckenberg |2|, Fuks [10], Boccara and Fuks [8], Belitsky and Ferrari [4], Belitsky, 
Krug, Neves and Schiitz [5J, Blank [HI [7], Helbing [14] . Wolfram [22] among many 
others. 

We introduce a slow-to-start traffic model which is continuous in space and time. 
Initially cars are distributed in the straight line M following a Poisson process of 
parameter A; cars may have speed 1 or 0. All cars start at zero speed and each car 
waits a random (delay) time exponentially distributed with mean 1 to to change its 
speed from to 1. Delay times of different cars are independent. After the delay time, 
each car moves at speed 1 until it collides with a stopped car to its left or forever if 
it does not collide. When car i collides with car i — 1, its speed drops immediately to 
zero and remains blocked until car i — 1 leaves. Then car i waits a further random time 
with exponential distribution before moving. And so on. At each time there are cars 
at speed 1 and cars at speed 0. We prove that if A < 1 every car will eventually be 
unblocked forever and the final relative positions of cars are distributed as a Poisson 
process of rate A. 

The main tool is a relation between the traffic model and a M/M/l stationary 
queuing system. Customers arrive at rate A according to a Poisson process to a single 
server whose service times are independent and exponentially distributed with mean 1. 
There exists a stationary version for the queue when A < 1. The stationary process 
is reversible so that its distribution is invariant by reversing the arrow of time. As a 
consequence the customer departure process is a Poisson process as the arrival process. 
This is the famous Burke's theorem. 

In Section [2] we define the process and give the main results. In Section [3] we 
construct a semi-infinite traffic model and define traffic cycles. In Section [4] we define 
the queue process and construct workload cycles. In Section [5] we relate traffic final 
relative car-positions and customer exit times in the queue. In Section [6] we use an 
approach of Thorisson to construct the stationary traffic process using cycles and in 
Section [7] we do it for the queue and conclude the proofs of the results. 
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2. Definition of model and main results 

We consider a system of cars moving in R from right to left. Cars have two possible 
velocities: either or 1. Cars are represented by points and cannot overpass. A car 
needs an exponential random time of mean 1 to pass from speed to speed 1. A typical 
car starts at speed 0, waits an exponentially distributed random time to change its speed 
to 1, travels at this velocity until it is blocked by another stopped car, waits the stopped 
car to leave, waits another exponential time to get again velocity 1 and so on. 

We construct a sequence (II, V) = ((II(t), V(t)); t > 0) of car trajectories, 
U(t) = (7Ti(t),i G Z) and V(t) = (vi(t),i G Z). For each i, 7T, : IR+ — > R is a piecewise 
linear function almost everywhere different iable; 7Tj(t) represents the position of car i at 
time t and Vi(t) G {0, 1} its velocity. The initial car positions are given by (yi, i G Z), 
with k 6 R and yi < yi + i for all i: set 7Tj(0) = yi. The initial velocities are all null: 
Vi = 0. The trajectories must satisfy the following properties: 

(i) 7Tj(0) = yt 

(ii) 7Tj(t) = —Vi(t) if Vi is continuous at t 

(iii) Vi(t) jumps from to 1 at rate 1 if iii(t) > 7Ti_i(t). 

(iv) Vi (t) = if m(t) =n-i(t) 

If instead of taking i in Z we consider a semi-infinite configuration of cars with 
initial positions yo < y\ < . . ., the trajectories can be constructed inductively. Car 
waits an exponential time of rate 1 and then assumes velocity 1; since it will not be 
blocked by any other car, it will continue at speed 1 for ever. Car 1 does the same, 
but if it collides with car (at yo), then it stops and after car leaves, it waits an 
extra exponential time and then goes, and so on. Let 7T, := (vrj(t), t > 0). If we think 
i as discrete time, the process of trajectories is Markovian in the sense that the law 
of the trajectory 7^ given 7Tj_i, . . . ,n depends only on 7Tj_i. Our first result says that 
it is possible to construct a spatially translation invariant traffic process if the initial 
positions of the cars are given by a Poisson process of rate A. 

Proposition 1. If (y,; % G Z) is a stationary Poisson process in M. with rate A and A < 1 
then there exists a spatially stationary version of the process IT = ((7Tj(i); t > 0), % G Z) 
with initial positions vrj(0) = yi for all i G Z. 

In the stationary Poisson process by convention yo < < y±, so that —yo, y\ and 
yi+i — yi for i ^ are iid random variables exponentially distributed. Let 

U := sup{s : Vi(s) = 0} 

(possibly equal to infinity) be the last time car i has velocity 0. From ti on, car i goes 
freely at speed 1. We say that car i is free at times t > ti. 

Proposition 2. Assume A < 1 and consider the traffic process U of Proposition^ For 
each i G Z, £, is finite almost surely. 



A slow-to-start traffic model related to a M/M/l queue 



4 



Let Di be the total delay of car i defined by 

A := \ - U - MU) - Vi)\ (3) 

Here —U is the displacement car i would have at time U if it started at speed 1 and had 
never been blocked and 7Tj(^) — y^ is the actual displacement by that time. The absolute 
value of the difference is the delay of car %. Call 

Si := yi + A (4) 

To understand the meaning of Sj, observe that if a car starts at position Sj at speed 
1 and it is never blocked, then its position at time t coincides with iTi(t) for all times 
* > U. 

Proposition 5. If A < 1 then {sf, % e Z} is a Poisson process of parameter A. 
Furthermore, 7Tj(i) — 7Tj_i(£) = s« — Sj_i /or t > maxjtj, 

The statement is that {s^; z e Z} as a subset of R is a Poisson process; the 
specification of the indexes may corrupt the Poisson property. These results say that if 
the initial position of cars are distributed according to a Poisson process and all with 
zero speed, then every car will eventually have velocity 1 and the relative car positions 
will be distributed according to a Poisson process. For this reason we call Sj the final 
relative position of car i. 

Let y G K and for initial position yi > y let r\ be the random time defined by 

n{y) : = sup((7r i )" 1 (t/)) 

this is last time the time car i is at position y. Let T(y) be the first time all cars crossing 
y after T(y) are free: 

T(y) := inf{rj(y); U < r m (y) for all m > i) 

Proposition 6. If X < 1 then T(y) < oo almost surely for all t/GR. 

3. Construction of trajectories for a semi-infinite initial car configuration 

We construct car trajectories for initial car positions = y < y\ < . . . and relate them 
to a M/M/l queue starting with an empty system at arrival of a customer. 
The trajectories will be defined as a function of a marked Poisson process 

(y,£) = ((s/i,6,m); *>o, o<m<i) (7) 

where y = (yi,i > 0) is a Poisson process on [0, oo) with rate A with a car added at 
the origin (that is, y — and (y*+i — yf,i > 0) are iid exponential random variables 
with mean 1/A) and the delay times £ = ((£j, m ; < m < i), i > 0) are iid exponential 
random variables with mean 1. The sequences £ and y are independent. These times 
are used as follows. Car i may collide with car i — 1 at sites yo, ■ ■ ■ , Vi-i', if the collision 
occurs at y m , then after car i — 1 leaves y m , car i waits £ i>m units of time before taking 
again speed 1. More rigorously, car starts at the origin and waits £ ,o units of time 
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to start moving at speed 1. Since there are no cars to the left of car 0, it will never be 
blocked and we define 

M / if < « < {„,„ 

The trajectory of car i is then defined as a function of the trajectory of car % — 1 and 
the waiting times < m < z) as follows. We shall define Aj im as the time car i 

arrives to y m and B^ m as the time car i departs from y m . Clearly A it7n < B im and if 
car i is not blocked at y m , then A iym = B iym . Let 

• 0, Bi i . j 

and then inductively assume Aj_ l m , Bi_ 1)m are defined for all < m < i — 1, as well as 
Ai )Tn and B iiTn . Then set 

Ai,m— i ■ B^ m -\- y m y m —\ (9) 

D I Ai m_i if -Aj.m-l > ^i-l.m-1 / 1n \ 

#i,m-l := < ' . ., 4 R UUj 

I -£>i-l,m-l + ?i,m-l 11 A i,m-l < ^i-l.m-l 

In words: the arrival of car i to site m — 1 occurs (y m — y m -i) time units after its 
departure from site m. The departure of car i from site y m -\ occurs immediately (at 
arrival time) if car % — 1 has already left or, £j im _i units of time after Sj_ 1>m _i, the 
departure of car z — 1 from site y m _i. 

The vector ((A iTn ; B iim ), < m < i) determines the trajectory (7Tj(t),t > 0): 

i 

id(t) :=^y k l{te (A\fc) B ijk )} + (y k + B ijk — t)l{te (B ifk , A,fc-i)} (H) 

fc=0 

where 1{-} is the indicator function of the set {•}. The total delay of car % defined in 
(J3]) satisfies 

i 

A = £(fli, fc - A ijfe ) (12) 

fc=0 

To stress the dependence on (y, £) we write A,fc(y,£), Bi ik (y,£), etc. Let s(y,£) = 
(soi s i) • • •)) where the final relative position Sj is defined as a function of and A as 
in (jlj). Let <x(y, £) = (ctq, cri, . . .) be the sequence defined by ctq := £ ,o and for i > 1 

{A,0 — if Sj-l > Hi, QgN 

otherwise 

Uj is called the final delay of car i. 

Lemma 14. cr = (cr , (Ti, . . .) is a sequence of iid random variables with exponential law 
of mean 1. Furthermore cr is independent of (y*; i > 0). 
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Proof Fix the trajectory of car i — There are two cases: either car i is blocked at by 
car % — 1 or not. In the first case B i0 = -Bi-i,o+£i,o> so that = £ ij0 which is independent 
of the trajectories (ir m ; < m < i — 1) and in particular of (o~ m ; < m < i — 1). In the 
second case the label of the leftmost blocking position of car % is given by 

K : = min{/c < i : Bi_ 1>k + f i)Jfc > B^ lfi - (y k - y )} 

here by convention B^i i = 0. K is a stopping time for (£i,i_ m ; m > 0); that is, the 
event {K = k} is a function of (£i,i- m ; < m < fc). But the dependence on ^ is only 
on the event {£ i>k > B^ lfi - (y k - y ) - Bi_ ljk }. Then, given this event, 

°~i = ii,K — — (yK — Vo) — Bi_i tK ) 

is exponentially distributed with mean one and independent of (n m ; < m < i — 1). 



Traffic Cycle Call 

X := mm{yi > : ^ > s^} (15) 

iV := mm{i > : y { > s^} (16) 

C := (((tt^O; * > 0), i e {0, . . . , N - 1}), TV, X) (17) 



We say that C is a cycle with length X and TV cars involved. The cycle C consists 
of a space interval X and TV car trajectories in the time interval [0, oo) with starting 
positions in [0, X); however, since car i is free after time B i0 , the trajectories are 
determined by the set ((n^t); t G [0,B ii0 ));i e {0, . . . , N — 1}) or, alternatively 
by the arrival/departure times of the TV cars to/from sites y ,...,yN-i given by 
((A,m, £i,m); o < i < TV - 1, o < m < i). 

The cycle C induces a stochastic process ^(C) = (Z y (C), y G [0,X)) given by the 
interval-valued vector (with dimension depending on y) 

Z V (C) := (irr\ y ) ; y < n < X ) (18) 

The fcth coordinate of Z y {C) contains the time interval spent at y by the fcth car to the 
right of y. Assuming this car has label % there are two cases: (a) if y = y m for some m, 
the interval is n^iy) = [Ai >m , Bi <m ] and (b) if y is not an initial car position the interval 
is a point, because car k will not be delayed at y. Z y (C) is a vector of zero length if 

y e [vn-i,x). 

The cycle C, its length X and its number o cars iV are functions of y = (y , yi, . . .) 
and £ = (f i>m ; % > 0, < m < i): 

C = C(y,0; X = X(y,|); N = N(y,£) 
Given a cycle C we can recover the length of the cycle, the number of cars involved, the 
initial car positions, the final delay times and the final relative car positions which are 
denoted 

N(C),X(C); y t (C), ^(C), Si (C), i = 0, . . . , N(C) - 1 . (19) 

The final relative car positions in the cycle coincide with the last passage through the 
origin: 

Si = B it0 , for % = 0, . . . , N - 1 (20) 




Figure 1. The first cycle has 4 cars. Car 1 collides with car and car 2 collides with 
car 1 at yo (the trajectories are drawn slightly separated at yo; the actual trajectories 
partially intersect at yo)- Car 3 does not collide with previous cars, but since its 
starting position 7/3 is to the right of S2, it still belongs to the cycle. Car 4 initial 
position 2/4 is to the right of S3 so that the new cycle starts at Si =2/4. 



4. FIFO M/M/l queue 

The FIFO M/M/l queue is a Markov process in N = {0, 1, . . .}. At rate A > customers 
arrive into the system, stay in line and one customer at a time is served at rate 1. This 
means that the service times are exponentially distributed with mean 1. The customers 
respect the time of arrival: first in, first out in the jargon of queuing theory. We use 
the following notation 

• fji arrival time of customer % 

• <7j service time of customer i. 

• §i exit time of customer i. 
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The queue size at time y can be constructed as a function of a marked Poisson 
process cr*); i G Z). The arrival process (y"j ; z G Z) is a Poisson process of 
parameter A. The service times (5$; i G Z) are mc2 random variables with exponential 
law of mean 1. The sequences are independent. 

The workload W y , y G R is the amount of service due by the server at time y. This 
is the time a customer arriving at time y has to wait to start to be served. This process 
is continuous to the right with limits to the left and piecewise linear. W y jumps at times 
yi, the arrival time of customer i by an amount <7j, the service time of this customer 
and it decreases continuously with derivative —1 until hitting 0, where stays until next 
arrival. The workload process must satisfy the following evolution equations: 

^y = -\{W y >0} fory^y,, zGZ (21) 
dy 

Wft=Wfc_ + ai (22) 

Construction of the workload process with an arrival to an empty system Let y = 
{Ho, Vii • • •) be a Poisson process of rate A as y, with y = and let <x = (a , &i, . . .) be 
a sequence of iid exponential random variables with mean 1 as er. Define 

W := a (23) 

and then, recursively, for i > 1: 

Wft := [Wy^ - (y t - y^)] + + cr, (24) 
W„ := [W^^ - (y - y G (y^x, y<) (25) 

To stress the dependence of (y, or) we denote W(y, 5") = (^(y, or); y > 0) the process 
defined by (T23|). (|24l) and (|25l). 

The exit time of customer i is defined by 

h := m + ^ (26) 

for 2 > 0. That is, the ith arrival time plus the workload of the server at arrival (included 
the ith service time). Define 

s(y, or) := (s , h, ■ ■ •) (27) 

Workload cycle Let X be the first time after time zero the workload jumps from 
to a positive value: 

X := inf{y > : W y - = 0, W y > 0} 

Define 

C:= (W y ;ye[0,X)) 
and let N be the number of arrivals during the cycle: 

N := max{z : y t G [0,X)} + 1 
(we add 1 to take account of the arrival at time 0). 
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Clearly C, X and N are function of (y, <x) : 

C = C(y,cr), X = X(y,cr), N = N(y,cr) 

Lemma 28. If X < 1 then X has finite expectation and there exists a stationary version 
of the workload process denoted W = (W y ; y G R). 

For a proof see Loynes [16J; we prove this later using cycles, following [20J. The 
proof of the following theorem can be found in Prabhu [IB] , theorem 8 (page 98) or 
Baccelli and Bremaud [JJ. 

Theorem 1 (Burke's Theorem). The exit times (sf, % G Z) of the stationary process W 
is a Poisson process with parameter A. 
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Figure 2. Workload evolution related with the traffic example of Figure 1. The cycle 
has 4 customers. The exit time S3 of customer 3 occurs before the time arrival j/4 of 
customer 4 so that a new cycle starts at S± = J/4- The service times &i coincide with 
the final delays Oi and the customer exit times Si coincide with the final relative car 
positions Si of Figure 1. 



5. Queue associated to traffic model 

The car trajectories generated by semi- infinite car positions y and car-delays £ defined in 
(J7J) generate final car-delays <x = <r(y, £) and last passages through the origin B i>0 (y, £). 
The queue generated by arrivals y and waiting times a produce a workload process 
W y (y, &). The relation between these processes is given by 

W y = (B ifi - y) + for y E [ Vi , y i+1 ) (29) 
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This follows from the definitions (j23j) - (|25j) and (fT3"|) . In the corresponding cycles this 
relation reads 

V \ if y^_i < y < y N { ' 

As a consequence of (1291) . the queue exit times s = s(y, cr) coincide with the final relative 
car-positions s. The length and number of elements in the respective cycles agree: 

Lemma 31. Let (y,£) be semi infinite car initial positions and delay times as in |?|). 
Then 

s(y,<r(y,£))=s(y,0 (32) 
X(y,*{y,£)) = X{y,£), N(y, <r(y,£)) = N(y,£) (33) 



Proof (1331) is a consequence of (1321) and the definitions. It suffices to prove (1321) for 
each cycle. a = £ ,o — Ah so that by (j3|, s = + cr = So by (126]) because W Q = oo- 
For 1 < i < N, 

Si = yi + W Vi = B ifi = Si 

by tfHD, d30D and (J2QD. 



6. Stationary traffic process 

Let ((y n , £ n )'i ^ G Z) be a iid sequence of Poisson processes and delay times with the 
same law as (y, £)• Let C n = C(y n ,$, n ) be the cycles generated by these variables. 
Then (C n ; n G Z) is a sequence of iid cycles with the same distribution as C . Let X n 
be the length of cycle C n , N n the number of cars involved with this cycle. By (1331 and 
Lemma [28] X n has finite expectation. Let 

{0 if n = fo if n = 

ELi X k if " > M n := J ^Li ^ if ^ > (34) 

Define the traffic process Z° = (Z°; y G R) by 

Z° := Z y _ Sn (C n ) for ye [^,5^), nGZ (35) 

as the process obtained by juxtaposing the cycles one after the other and putting the 
beginning of cycle at the origin; recall the definition of Z y (C) in (118j) . In the process 
Z° the cars of a cycle are not blocked by cars of the previous cycle; in particular, all 
cars starting at positions S n , n G Z, will be free after waiting an exponential time of 
mean 1. 

Let U be a uniform random variable in [0, 1] independent of Z°. Let P° be the law 
of (((y n , £ n )i n ^ ^)>^0 an d E° the corresponding expectation. As Z° is a function 
°f ((yn,^n)s n e ^)) we can a l so think that P° is the law of (Z°,U). Let # r be the 
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translation operator defined by (8 r Z) y = Z y _ r . Define Z as a function of Z° and U as 
the process Z° "as seen from —(1 — U)X ": 

Z := O-a-wxgZ yeR (36) 
and define the law P by size biasing cycle zero: 

dF := J^dF" (37) 

E°X v ; 



The following result is in Theorem 4.1 of Chapter 8 of Thorisson [20]; see also Figure 
2.1 in that chapter for a better understanding of the relation between Z and Z°. 

Lemma 38. The law of Z under P is stationary. 

For test functions / (from the space where Z is defined to M.) the law of Z under P 
satisfies 



Ef(Z) 



E°[X /(^ Xo(1 _ c/) Z 



E°X 

where Xq is the lenght of cycle of Z°. To obtain a sample of Z under P, first sample 
a process with a cycle starting at the origin with the size biased law P, then translate 
the origin to a point uniformly distributed in cycle 0. 

Since the traffic process Z is constructed as juxtaposition of cycles, the initial car 
positions, the final delays and the final relative car positions can be recovered as follows 
using the notation of ffT9l) : 

S = sup{y < : Z y _ = 0, Z y > 0}, (39) 

S n = inf {y > 5 n _! : Z y _ = 0, Z y > 0}, n > 0, (40) 

S n = sup{y < S n+1 : Z y _ = 0, Z y > 0}, n < ; (41) 

X n = S n+1 - S n , C n = (Z y ^ Sn ; y e [S n , S n + X n )), N n = N(C n ) ; (42) 

Vk = S n + yk-M n (Cn), Cfc = 0~k-M n (C n ) , Sfc = S n + Sk-M n (Cn), (43) 

for M n <k< M n+1 . (44) 

Denote y(Z), cr(Z) and s(Z) the initial car-positions, the final delays and final relative 
car-positions of the process Z. 

Lemma 45. The law of y(Z) and tr(Z) under P are stationary and y{Z) and cr{Z) 
are independent. In particular y(Z) is a Poisson process of parameter A and cr(Z) is a 
sequence of iid exponential random variables of mean 1. Furthermore s(Z) is a Poisson 
process of parameter X. 

Stationarity of y(Z) and cr(Z) is immediate consequence of stationarity of Z. The 
other properties follow from the stationary construction of the queue and relation (J3"2l 
as shown in the next section. 
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Trajectory extension The trajectory of the car starting at y^ G [S n , S n+ i) is defined 
only for t G [0, -B ij0 (C n )]; see ffTB"]) . We extend the trajectory by just continuing at speed 
1 from this point on: 

n . {t) = | S n + 7rl Mn (t) if Vl G [S n , S n+1 ), t G [0, B i>0 (C n )) 

) S n — t + Bi-M n ,o(C n ) if ?/t G [S n , S n +i), t > B i: o(C n ) 

where 7r™ is the trajectory of the ith particle of cycle C n . Using definitions (146]) and (jTSI) 
the corresponding process Z has a countable number of coordinates at each position y; 
the kth coordinate indicates the interval of time spent at y by the kth car initially to 
the right of y. We abuse notation and continue calling Z this process. Notice that y(Z), 
cr(Z) and s(Z) remain unchanged and Lemma 1451 holds for this extension. 

7. Stationary workload process 

Assume A < 1 and let (y n ,£ n ) be the sequence introduced in the beginning of 
previous section. Let cr n = <x(y n ,£ n ) as defined in ( [131) . By ( 1331) the workload cycles 
C n = C(y n , cr n ) have the same length as the traffic cycles C n : X n = X n and the number 
of cars in cycle C n is the same as the number of customers in cycle C n : N n = N n . Let 
S n as in flM]) and define W° = (W°; y G R) by 

W° = W y _ Sn (Cn) forye[S° n ,S° n+1 ), neZ (47) 

and the process W = (W y ; y G R) by 

W y = W° H1 _ U)XS yeR (48) 
where U is the same variable used in (|36|) . As before, for P given by (1371) . 
Lemma 49. Under ¥ the law ofW is stationary. 

W is a stationary M/M/l queue with arrivals y(W) and departures §(1^). Hence 
the arrival process y(W) is a stationary Poisson process of rate A in R. By Burke's 
Theorem (TJ the same is true for the departure process s(VK). 

The stationary workload process W and the stationary traffic process Z are 
constructed in the same space (as function of ((y„,£ n ); n G Z),Z7) so that they have 
exactly the same cycles and the initial and relative final car positions of Z coincide 
respectively with the arrival and departure process of W: 

y(Z)=y(W), s(Z)=~s(W) (50) 

This finishes the proof of Lemma [45] and Proposition [5j 

8. Final remarks 



When the initial density of cars is smaller than the inverse of the delay time, A < 1 
the process is called subcritical. We have described how a stationary configuration 
of initial car positions organize the departure from speed zero to speed 1 under the 
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rule slow-to-start in the subcritical case. The method relates the space-stationary one- 
dimensional slow-to-start traffic model with the workload process of a M/M/l time- 
stationary queuing system. If the initial position of cars is Poisson then the final relative 
position of free cars is also Poisson. The same method shows that if the initial position 
of cars y is a stationary ergodic process with density A < 1 (density is the mean number 
of cars per unit length) then the final relative car position process coincides with the 
departure process of the queue with arrival process y. 

The supercritical case merits to be investigated. When A > 1 cars form "traffic 
jams" on a subset of the initial positions y, depending on time. As time grows, the 
density of this set goes to zero and the number of cars per traffic jam increases. We 
agree with an anonymous referee who conjectures that as t goes to infinity the set of 
traffic jams suitably rescaled converges to a Poisson process of rate 1. The critical case 
A = 1 should also show traffic jams. 

Another model to research might include spontaneous stops at some rate v keeping 
the slow-to-start rule. These kind of models will be closer to those studied by Nagel 
and Schreckenberg and Gray and Griffeath. 

We are investigating the same phenomena for cellular automata in Z. The approach 
can be applied but there are complications coming from the hard core interaction. The 
results are similar. 
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